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Magnetic fields during high redshift structure formation 
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We explore the amplification of magnetic fields in the high-redshift Universe. For this purpose, we perform high-resolution 
cosmological simulations following the formation of primordial halos with ~ 10 7 M©, revealing the presence of turbu- 
lent structures and complex morphologies at resolutions of at least 32 cells per Jeans length. Employing a turbulence 
subgrid- scale model, we quantify the amount of unresolved turbulence and show that the resulting turbulent viscosity has 
a significant impact on the gas morphology, suppressing the formation of low-mass clumps. We further demonstrate that 
such turbulence implies the efficient amplification of magnetic fields via the small-scale dynamo. We discuss the proper- 
ties of the dynamo in the kinematic and non-linear regime, and explore the resulting magnetic field amplification during 
primordial star formation. We show that field strengths of ~ 10 -5 G can be expected at number densities of ^ 5 cm -3 . 
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1 Introduction 



In the present-day Universe, m agnetic fields ha ve been de- 
tected in galaxies jBeckL[2QQ4h . galax y clusters (Kim et al 



Peters et all 20 


L2l: Schleicher et all 201 


0l Schober et al.. 


2012a: ISur et al. 


liom 2010: Turk etal.. 


2012b. 



A central condition for this scenario is th e presence 



1990h and local dwarf galaxi es ( Chyzy et all 2011 



Hees en et all l201ll: iKepley etall l201ll) . Evenfor the i 



of turb ulence. The latter has been reported by Abel et al 
d2002l) for the first star-formi ng minihalos, and iGreif et al. 



tergalactic medium, more speculative o bservations sug 



gest t he presence of magnet i c fields ( Broderick et al 
20121: iMiniati & Elvrvl l2012l: TNeronov & VovkL 12010 : 



l2008h and lWise & Abel d2008l) for larger primordial halos. 
However, such massive halos were only exp lored with a res 
olution of ~ 16 cells per Jeans length, while iFederrath et al 



Takahashiet alll2012l:lTavecchio etaD.l201ll:lYuksel et al 
20121) . Magnetic fields have further been inferred in high- 
redshift galaxies via the observed correlation of quasar ro- 
tation measures with t he number of galaxi e s along the line- 



of- sig ht until z ~ 2 (IBernet et all 120081: iKronberg et al 



2008), as well as through the observed far-infr ared - ra- 
dio co rrelation, which is established until z ~ 4 (IMurphy . 
2009h . 



We propose here that strong magnetic fields are pro- 
duced as a result of turbulent magnetic field amplifica- 
tion already at high redshift. The small-scale dynamo 
efficiently converts turbulent into magnetic energy, on 
timescales considerably small e r than t he dynamic a l time 
(iBrandenburg & Subramaniarl 120051: iKazantsevl Il968 : 
Schober et all 2012bh . As a result, the dynamo process will 



likely produce strong magnetic fi elds during the forma- 
tion of the first stars and galaxies (iFederrath et all bOllbt 
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have shown that at least 32 cells per Jeans length 
are required in order to obtain converged turbulent energies. 

In this paper, we present a detailed study of the turbulent 
properties of massive primordial halos, exploring the ef- 
fects of resolved turbulence by systematically increasing the 
numerical resolution per Jeans length, as well as the unre- 
so lved turbulence using the su bgrid- scale turbulence model 
by ISchmidt & Federrathl d201 ll) . We then discuss the ampli- 
fication of magnetic fields in the kinematic and non-linear 
regime. Our results indicate the likely presence of strong 
magnetic fields in the high-redshift Universe. 

2 Turbulence during the formation of 
protogalaxies 

To explore the generation of turbulence during the for- 
mation of protogalaxies, we perform cosmological hydro- 
dynamics simulations with the adaptive -mesh refinement 
(AMR) code EnzcQ dO'Shea etall 12004b . We explore the 
evolution of a cosmological box of 1 Mpc/h with a root grid 
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resolution of 128 3 and two initial nested 128 3 grids cen- 
tered on the most mass ive halo, each inc reasing the reso- 
lution by a factor of 2 dLatif et all I2012|) . We include pri- 
mordial chemistry following the non-equilibrium evolution 
of H, H + , e~, H~, and H2, along with a photodisso- 
ciating background corresponding to J21 = 1000, where 
J21 is given in units of erg s _1 cm -2 Hz -1 sr _1 . Our ver- 
sion of the code i ncludes a subgrid-scale model for unre- 
solv ed turbulence (Ma ier et aD.l2009l;ISchmidt & Federrathl 
12011b . accounting for turbulent pressure, diffusion and in- 
cluding an explicit modeling of the turbulent energy dissi- 
pation. Our main refinement criterion is the number of cells 
per Jeans length, where we explore values of 16, 32 and 64. 
We stop our simulations when the highest refinement level 
27 is reached. 

The halo collapses at a redshift of 11.9 with a total mass 
of 8.06 x 10 6 M . The resulting density distributions in 
the central 500 AU are given in Fig. [T] for different res- 
olutions per Jeans length, comparing the results with and 
without the subgrid-scale model. While the low-resolution 
runs may give the impression of rather smooth density struc- 
tures, a transition to turbulent, complex structures becomes 
visible for a resolution of 32 cells per Jeans length, which 
is more pronounced when 64 cells are adopted. This trend 
occurs both in the runs with and without the subgrid-scale 
model. The morphologies may however change consider- 
ably in the presence of the subgrid model, especially for 
the high-resolution results. The latter provides an indication 
that turbulent structures are not converged even in the high- 
est resolution runs. For the radial profiles of the halo, we do 
however find convergence, as illustrated in Fig. 

We have examined the trends described here for a set 
of three massive halos, as described in more detail by 



Latif et al.l (|2012|) . For these halos, we calculated the prop- 



erties of gas clu mps within the cen tral 500 AU using the 
clump finder of Wil liams et al] (|1994|) . As shown in Fig. [3 
the latter exhibit a rather similar power-law behavior in all 
halos considered. We further find that the formation of low- 
mass clumps is suppressed in the presence of the subgrid- 
scale turbulence model. 



3 Turbulent amplification of magnetic fields 

The turbulence present in these halos may efficiently am- 
plifiy even weak magnetic fields via the small-scale dy- 
nam o. Appropriate see ds may result from the Bierm ann bat- 
tery (Biermann, 1 195(1 the W eibel instability ( Lazar et al 



2009; Medvedev et all 120041) . thermal plasma fluctuations 



dSchlickeiserl l2012|) or even the p r e-recombination Uni 



verse dBanerjee & Jedamziki 120031: iGrasso & Rubinstein! 



2001) 



The small-scale dynamo has been explored 
with numerical simulations and analytical models 
( Brandenburg & Subramaniani |2005|) . An analyti- 
cal trea tment is possible in terms of the Ka z antsey 
model dBoldvrev & Cattaneol 12004 iKazantsevl 1 19681: 
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Fig. 2 The radial profiles of our halo at z = 1 1 .9 (top left: 
gas density, top right: gas temperature, bottom left: total en- 
ergy, bottom right: subgrid-scale energy) for different reso- 
lutions per Jeans length, with and without the subgrid-scale 
model. Unlike the detailed morphology, we find that the ra- 
dial profiles already converge for moderate resolutions per 
Jeans length dLatif et al Iko 12h . 



Schober et all l2012bl: ISubramaniani Il998l) . describing 
the dynamo amplification as the advection of passive 
vectors in a velocity field consisting of homogeneous and 
isotropic turbulence. The model describes magnetic field 
amplification in the kinematic regime, before backreactions 
occur. We note that two phases are relevant for the origin 
of high-redshift magnetic fields: While the kinematic 
phase leads to a fast exponential growth on timescales 
comparable to the eddy-turnover time in the viscous range, 
the non-linear growth leads to the formation of larger- scale 
magnetic fields. The overall process is expected to saturate 
within a few eddy-turnover times. 

3.1 The kinematic regime 

We start our considerations with the kinematic regime of 
dynamo amplification. The latter is described in terms of 
the induction equation, which is given as 



dB 

— = VxvxB-r?VxVxB, 
dt 1 



(1) 



where B is the magnetic field, v is the fluid velocity 
and r\ the magnetic diffusivity. In the absence of helicity, 
the correlation function of the magnetic fields, M^ (r) = 
(Bi(0)Bj(r)), can be decomposed into a transversal and 
longitudinal component: 



Mij(r) = fa - -2 



'* r ^M N (r) + r 4fM L (r). (2) 



As a result of the constraint equation, V • B = 0, it can be 
shown that 



MN = i§- r ^ ML{r) ^ 



(3) 
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Fig. 1 The density distribution in the central 500 AU with and without the turbulent subgrid model, and for different 
resolutions per Jeans length (top to bottom: 16, 32 and 64 cells per Jeans length, left: standard hydro, right: hydro+subgrid- 
scale mode). We find that increasing the resolution per Jeans length considerably changes the morphology of the halo, as 
we start resolving turbulent structures at high resol ution. The subgrid -scale model favors the formation of more compact 
structures and reduces the number of density peaks dLatif et all 1201 2|) . 
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Model and reference 



label 



T (Pm -> oo) 



Kolmogorov (1941) 

Intermittency of Kolmogorov turbulence (She & Leveaue, 1994) 
Driven supersonic MHD-turbulence (Boldvrev et al., 2002) 
Observation in molecular clouds (Larson . 1981) 
Solenoidal forcing of the turbulence (Federrath et al., 2010) 
Compressive forcing of the turbulence (Federrath et al., 2010) 
Observations in molecular clouds (Ossenkopf & Mac Low, 2002) 
Burgers (1948) 



K41 
SL94 
BNP02 
L81 

FRKSM10 
FRKSM10 
OM02 
B48 



1/3 
0.35 
0.37 
0.38 
0.43 
0.47 
0.47 
1/2 



31 Re 1 / 2 

36 

0.94 Re 48 
0.84 Re 46 
0.79 Re 45 
0.54 Re 40 
0.34 Re 36 
0.34 Re 36 
^ Re 1/3 

60 



Table 1 The nor malized growth rate of the small-scale dynamo T in the limit of infinite magnetic Prandtl number 
( Schober et al 112012b). We show our results for different types of turbulence, which are characterized by the exponent $ 
of the slope of the turbulent velocity spectrum, v(£) oc . The extreme values of $ are 1/3 for Kolmogorov turbulence and 
1 /2 for Burgers turbulence. 
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Fig. 3 The properties of gas clumps in the central 500 AU 
for three different halos, calculated with and without the 
subgrid- scale model and for a resolution of 64 cells per 
Jeans length. At this stage, the clumps are still unbound, 
and many of their properties follow characteristic power- 
laws. We note that the number of low-mass clumps is 
strongly reduced in the presence of the subgrid- scale model 



(ILatif et all l2012|) . Colors correspond to different halos, 



triangles denote simulations with, diamonds without the 
subgrid- scale model. 



The same decomposition can be performed for the turbu- 
lence corr elation function Tg (r) = (vi(0)vj(r)). As de- 
scribed by ISchober et alJ d2012bl) . we adopt the following 
parametrizations for the longitudinal correlation function, 

'l_ Re (i-»/(i+*)(£) 2 ) <r<4 

t v < r < L 



ar 1 ) 



L < r, 



(4) 



and the transversal correlation function 



T N (r) = 



/ r 10+l' 



(f) 2 ) 







< r < t v 

i v <r < L 

L < r, 

(5) 

with c = VL/3 and U = (21 - 38tf)/5. The turbulent 
slope is given as d = 1 /3 for Kolmogorov and d = 1 /2 for 
Bur gers. The results for a larger set of models were derived 
by lSchober et al.1 (12012b|) and are given in Table [T] 
The latter are obtained employing an ansatz 

M L (r,t) = ^(r)e 2r V(r 2 V^(r)), (6) 

where t denotes the time coordinate, T can be interpreted as 
the growth rate of the correlation function, and n{r) a func- 
tion describing the turbulent diffusion. From this ansatz, it 
is possible to obtain the so-called Kazantsev equation, 



dr 2 



+ U(r)V(r) = -r*(r), (7) 



which has the same form as the time-independent quantum- 
mechanical Schrodinger equation with a potential U(r) 
given as 

(^) 2 



U(r) 



2Tn 



2(T L -T n -k) 



(8) 



It can thus be solved using quantum-mechanical meth- 
ods, in particular the so-called WKB- approximation. Intro- 
ducing the normalized growth rate f as 

f = (9) 



© 2013 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



www.an-journal.org 



Astron. Nachr. / AN (2013) 



5 



if) - 



1 r 



K41: tf=l/3 ■ 

SL94: $ = pig 
BNPC2: i? - 0.37 

J* LSI: = 0.38 

f > . - ■ FItXSM LO(syl)- = O.^J 

. - " ' FRKSM10 (comp ) ,OM.U2 : = 0.47 ■ ■«■ ■«■ ■ 

B4S: 1/2 ■ 

ml i M I id - i i mill i mid i iimiJ i ■imiJ 



10- 1 10° 10 1 10 2 10 3 Itjf 10* 

Re 

Fig. 4 The normalized growth rates of the small-scale dy- 
namo for different tur bulence models a s a function of the 
Reynols number dSchober et all 12012b ). The labels are de- 
fined in Table [U 



Schober et al.l (I2012b|) have derived the resulting ampli- 
fication rates in the limit of high magnetic Prandtl numbers 
Pm ^> 1, defined as the ratio of kinetic and magnetic diffu- 
sivity. This limiting case is indeed appropriate for the in- 
terstellar and intergalactic medium, while stars and plan- 
ets typically have magnetic Prandtl numbers of Pm <C 1. 
The resulting normalized growth rates are given in Table [T] 
The normalized growth rate is plotted as a function of the 
Reynolds number in Fig. H We find that the type of turbu- 
lence, and in particular its spectral slope has a major influ- 
ence on the resulting amplification rate, which can be re- 
duced by a few orders of magnitude. Nevertheless, we note 
that these growth rates are still orders of magnitudes smaller 
than the dynamical timescale of the system, implying that 
magnetic field amplification may still pro ceed in a highly 



efficient fashion in any of these cases (e.g. iFederrath et al 
201 id) . 



3.2 The non-linear regime 

After saturation has occured on the viscous scale, where am- 
plification originally is fastest, the magnetic field will be 
strong enough to prev ent further growth on that s cale. How- 
ever, as described by ISchekochihin et al.1 (I2002|) . magnetic 
field amplification will continue on larger scales £, with the 
typical eddy-turnover time V£ oc I® on that scale. For Kol- 
mogorov turbulence, the latter leads to a phase of linear 
growth of the magnetic energy, and one can show that Burg- 
ers turbulence leads to a quadratic increase of magnetic en- 
ergy with time. During that phase, the length scale of the 
peak magnetic fluctuations increases considerably, from the 
viscous scale up to the injection scale of turbulence. At that 
point, the system is saturated. At low to moderate Mach 
numbers, the magnetic energy consists of 30 — 60% of the 
turbulent energy, while at high Mac h numbers, smaller satu - 
ration values of ~ 5% are possible ( Federrath et all Eo 1 1 a|) . 



The time estimates required to reach saturation in 
the non-linear regime range from a few to a few dozen 
eddy- turnover times dBeresnyakl 1201 2l: ISchekochihin et al ' 



|2002|) . As a result, strong magnetic fields can be expected 
shortly after the formation of the first virialized objects due 
to turbulent amplification. 

3.3 Magnetic fields during primordial star formation 

In order to explore the consequences of turbulent mag- 
netic field amplifica tion for primordial star formation, 
ISchober et al.l d2012a|) have implemented a magnetic field 
amplification model within the one-zone framework fo r 
primordial collapse developed by iGlover & Savinl d2009|) . 
Their model assumes gravitational collapse on the free-fall 
timescale, solves the rate equations for primordial species 
along with the cooling/heating balance for the gas temper- 
ature (see Fig. \5\ for the thermal evolution). For the ini- 
tial magnetic field, we adopted a typical value of 10 -20 G, 
corresponding to the magnetic field strength from a Bier- 
mann battery. For the initial growth of the field during 
the kine matic regime, we em ploy the growth rates de- 
rived by Schob er et al. I d2012bh . In this calculation, we as- 
sumed that turbulence is injected on the Jeans s cale as 
a result of gravitational coll apse dFederrath et all 1201 ltJ: 
iKlessen & Hennebellell20ici|) . with typical turbulent veloc- 
ities corresponding to the typical sound speed of the gas. 
After saturation occurs on the viscous scale, the evolution 
of the magnetic field is ca lculated employing the model of 
Schekochihin et al.l ( 2002 ). We assume that the magnetic 



field is saturated when equipartition with turbulent energy 
is reached. 

We find that the magnetic field strength rises quickly 
during the kinematic regime, as the characteristic time 
scales are considerably smaller than the dynamical time. In 
the non-linear regime, the growth rate decreases when the 
saturation level is approached, which occurs after only a 
modest increase in density. The resulting field strength at 
densities of a few atoms per cm 3 is thus of the order of 
10 -5 G. The evolution of the field strength both for Kol- 
mogorov and Burgers type turbulence is illustrated in Fig.O 

Overall, we conclude that magnetic field amplification 
is highly efficient in the early Universe, resulting from the 
ubiquity of turbulence and the efficient amplification via 
the small-scale dynamo. We expect that such turbulent field 
structures can be probed in high-redshift starbursts with the 
Square Kilometre Arra}@. 

Acknowledgements. DRGS, JS, RSK and SB acknowledge fund- 
ing from the Deutsche Forschungsgemeinschaft (DFG) in the 
Schwerpunktprogramm SPP 1573 "Physics of the Interstellar 
Medium" under grants KL 1358/14-1 and SCHL 1964/1-1. DRGS, 
JN and WS thank for funding via the SFB 963/1 on "Astrophysi- 
cal flow instabilities and turbulence". JS acknowledges the support 
by IMPRS HD, the HGSFP and the SFB 881 "The Milky Way 
System". CF thanks for an ARC Discovery Projects Fellowship 
(DP110102191). 



SKA website: http://www.ska.ac.za/ 



www.an-journal.org 



© 2013 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 



6 



Schleicher et al.: Magnetic fields at high redshift 



References 




Fig. 5 The thermal evolution in primor dial gas dur- 
ing gravitational collapse, as calc ulated by ISchober et al 



(2012a) with the one-zone model of lGlover & Savin d2QQ9b . 
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